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Abstract 

We derive closed recursion equations for the symmetric polynomials occuring 
in the form factors of D^^^ affine Toda field theories. These equations follow 
from kinematical- and bound state residue equations for the full form factor. 
We also discuss the equations arising from second and third order forward 
channel poles of the S-matrix. The highly symmetric case of D'^^ form factors 
is treated in detail. We calculate explicitly cases with a few particles involved. 
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1 Introduction 



Form factors are matrix elements of local operators in a quantum field theory. Knowing 
them explicitly or at least perturbatively gives a deep knowledge of the quantum structure 
of a given classical Lagrangian field theory. For instance it might then in principle be 
possible to classify the local operator content of the theory or to calculate correlation 
functions. However, for a generic field theory it is a difficult problem to calculate these 
form factors. 

In the class of two-dimensional field theories one is in a better situation. It is known 
that there exist models which enjoy the property that scattering is factorizable, i.e. any 
n-particle scattering process can be decomposed into a product of two-particle scattering 
processes. 

For theories whose scattering is factorizable it can be shown that the form factors are 
subject to a set of axioms [|l|. The axioms (equations) provide a machinery which in 
principle enables one to calculate form factors explicitly. 

Within the last years based mainly on the studies |1| form factors have been studied for 
some diagonal scattering theories |0, |, ||, ^, H, ^, [l^ . Among them a complete solution 
to the form factor equations was obtained for the scaling Lee- Yang- model (minimal ^2^^) 
IP and for the sinh-Gordon model 0. However, these models are simple in the sense 
that they do involve only one type of particle or that the S-matrices of the corresponding 
theories do only have poles of first order. We note that only these first order poles are 
covered directly by the axioms mentioned above |Q. 

Recently the two particle form factors for the magnetic perturbation of the Ising model 
have been investigated This work is the first which treats a model with several 



species of particles and an S-matrix with higher order poles. 

Moreover some form factors in A^^ affine Toda field theories were studied in [1121. These 



theories do generically have more than one species of particles and the S-matrices have 
poles up to second order. Based on the analysis in [0 we would like to address the 
problem of form factors in D!^'^ affine Toda field theories. 

In order to make this paper a bit more self contained let us define affine Toda field theories 
(ATFT). They are given at the classical level by the following Lagrangian. 

1 

C = -<9^$ ■ 9"$ - V($), V($) = —J2n^exp{/3a^ ■ (1) 
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# = ($1, $2, • • • , ^r) is a vector of r scalar fields, and rii denote for i = 1,2, ... ,r 
the simple root vectors and the Coxeter labels of the corresponding root of a simple Lie 
algebra We consider only Lie algebras of type A, D, and E. 



The real constant m sets the mass scale of the theory while /3, which is assumed to be 
real throughout this paper, is the coupling of the theory. In fact it turns out that as long 
as we do not consider the perturbative structure we can work for our purposes with an 
"effective coupling" []T1| 



The exact S-matrices for ADE-ATFTs have been given in ||I5[ and serve as the main input 
for the form factor equations. It is known that S-matrices for the A-series possess poles 
of at most order two, while for the D-series they do have poles up to order four. In the 
latter case poles of first and third order correspond to forward channel physical particles. 
We will come to that point below. A general reference on this might be [|16], [1^. The 



form factor problem for these theories has first been addressed in 
In this paper we are going to shed some light on the form factors for D^p-ATFT. As it was 
already mentioned they provide an interesting complication of the studies quoted above 
since they do contain several different particle species and have S-matrices which contain 
higher order poles. 

The outline of the paper is the following. In section 2 we define form factors and review 



the axioms which they have to obey and some results of |T2[ which will be needed for our 
study. Some well known facts about S-matrices and fusings in D!^^ are recalled in section 
3. We then give a closed form of the kinematical and bound state residue equations and 
discuss the influence of the higher order poles to the form factor equations. In section 
4 the case of D'^^ will be discussed in detail. We point out how the symmetries of the 
Dynkin diagram enter the form factor equations and analyse in detail the poles of second 
and third order which are present in this particular model. We are going to calculate some 
examples of low particle form factors in D^^^ in section 5. The last section is devoted to 
comments and discussions. 
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2 Some facts about the form factor bootstrap 



As mentioned in the introduction a form factor is a matrix element corresponding to a 
local operator O. We take this matrix element in a special form which is commonly used 
in the literature: 

F,„„aM, ■ ■ ■ ,0n) ■■= {0\OmOu ■ ■ ■ ,0n). (3) 

In this expression 9i labels the rapidity of the particle of species a,, and |0) denotes the 
physical vacuum of the theory. The fact that we take the matrix element of O at the 
origin is simply a matter of convenience. 

The form factors are known to be subject to four axioms The first two of them are 
the so called Watson's equations. 

P'ai...aiai^i...an (^1; • • • ; ^i) ^j+l; • • • i ^n) Sa^ai^i i^i ^i+l)-fai ...ai+iai...a„ (^1; • • • ) ^i+l) ^i) • • • ) ^n); 

Fa^a2...a„{dl + '^'^h92,...,9n) = -^aa ...a„ai (6*2 , • • • , 6'n, 6*1 ) • 

(4) 

Needless to mention that Sab{0) denotes the S-matrix element. 

Due to the structure of the theory the form factors do contain poles. There are two kinds 
of them. The first one is due to the kinematics. The following equation can be explained 
by a field theoretical reasoning [0 and links a particle of species a with its antiparticle a. 

— i rese/=e+j^ Faadi...d,X^' ■, 6', 6^1, ... , On) = ■ ■ ■ ,0n) I 1 — H '^"'^j ~ ^i) j • (^) 

We will refer to this equation as kinematical residue equation. 

The second type of singularity is due to the possible bound states of the theory. As long as 
the singularities of the S-matrix which lead to bound states are of first order the following 
bound state residue equation for the fusion a + b —>■ c ai fusion angle 6'^^ and with on shell 
three-point vertex 

irj' = -I Te^e=^el,Sab{e)^ (6) 

holds 

- i TeS0'=0+igc^ Fabdi...d,Xd', OjOl,..., On) = ^'abFcdi...d„iO + iObcJ ^1, ■ ■ ■ , On). (7) 
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We have set ^ = tt — ^. 

The S-matrix of D^^-'-affine Toda theories is known to possess poles of order 2, 3 and 4 as 
well. The residue equations corresponding to these higher order poles need to be clarified 
and will be discussed below. 

The form factors are in general meromorphic functions in n variables. It is possible to 
split the full factor into a minimal part which is analytic in the strip < Im^^ < 27r 
and has no zeros in < Im^ < 27r. It can be shown in generalization of a theorem in 
that a solution of Watson's equations takes the form 

Fa,...aM, ■■■On)= K,....aM, ■ ■ ■ On) W F^^ie. - 6,), (8) 

i<j 

where it has to hold that F^^^^iO) = Sa,a,{e)F^^{-e) and F^^iO + 2m) = F^^^^-O). 
These minimal form factors are well known for ADE-ATFT. They can be presented in an 



integral form or equivalently in an infinite product expansion of F-functions (see e.g. 

my 

Now that these are known the work to be done is to solve for the form factor equations 
for the object K which contains poles and zeros in the above mentioned region. By 
the construction K automatically satisfies (^) by requiring some obvious monodromy 
properties. 

It has been shown in [|12| that from the kinematical (^) and bound state (|^ residue 
equations we get the following equations for K. 

—i rese'=9+j^ Kaadi...d„{0' , 9,9i, . . . , 9^) = 



Here we have set 



U{e)-' = n (10) 

xeAab 



This notation is taken from and will be explained in the beginning of the next section 



(0). Here we would only like to state that the S-matrix is given by 

Sa,{e) = (11) 

This form should be compared with the expressions for the S-matrix at the beginning of 
the next section. 

The bound state residue equation is 
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—i TeSg'=g+igc^ -K'a6di...d„ (6*', 6', 6^1, ■ ■ ■ , 6'„) — , , 

= r^,i^.-.,...„(^"+ ^1, ■ ■ ■ , On) n "=1 Kt;d,io + - o,)/Frm,)- ^ ' 

The object A is via Watson's equation very closely related to the S-matrix bootstrap 
equation. Since its explicit form will be of importance later we mention: 

= F^^Wm " -^J- \'^ac-x}+(e) 11 {x-u^^)+^e)- (13) 

■^—ac be 

The ^j^'s (^^^ = TT — ^^5^) are the fusion angles of the theory and we set u^^ = with 
h being the Coxeter number. 

We will give details on this notation at the beginning of the next section. 

We have now reduced the problem of calculating form factors in ATFT to that of making 
a proper ansatz for K and solving the two equations (|^) and ([l^). In cases where the 



S-matrix has higher order poles we do have to supplement these equations with additional 
ones. 



3 Form Factors in D^^^ Toda theory 



Let us recall some facts about D^^^^^-ATFT These theories contain n kinds of particles 
which will be labeled by the elements of the set {1,2, . . . ,n — 2, s, s'}. While the particles 
labeled by ordinary numbers in this set correspond to the points on the straight line, s 
and s' correspond to the particle upside and downside of the fishtail part of the Dynkin 
diagram of Dn- 

In this class of ATFT we have to distinguish between n odd and n even. In the even case 
all particles are self-conjugate while in the odd case we have s = s'. 
Let h = 2{n — 1) denote the Coxeter number of the theory in question. The masses of 
the particles are then given by 

ml = ml = 2171"^, = 8m^ sin^(^), k = 1,2, . . . ,n - 2. (14) 
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sets the mass scale of the theory. We remark that the case of d'^I'* is special in the 
sense that here we do have three particles {1, s, s'} of mass A/2m and one heavy particle 
{2} of mass \/Qm which corresponds to the central point of the Dynkin diagram of D4. 

The exact S-matrices and the fusion angles for the D^^^-theories have been given in |T5 



We recall them for convenience. In doing that we borrow the following notations from 
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_,„,:= Isinhiffl + ^r). (ry,y.= ^J±^, (15) 

and 



HQ) 



M (r + l)+(e)(r-l)+(g) _ 

V )+W - J——. m 7— m ' V )W - -rV\ ■ 

(r + 1 - i^j+(0)(r - 1 + {-r)+{e) 

Using these notations the S-matrices for n even are given by 

a+6-1 2a-2 

Sab{e)= W {p)ie){h-p)^e)^ Ssa{e) = Ss'a{e)= X{{n-a + p)^e). (17) 

\a-b\+l 
stGp2 stcp2 

with a, 6 G 1, 2, . . . , n — 2, and 

Sssie) = s,,A9) = n (p) w, SsAO) = U (p) w (is) 

1 3 
Step4 stcp4 



For odd n similar expressions are obtained [|T5] 



We denote by Aab the set of integers actually appearing in the brackets of the correspond- 
ing S- matrix above and by Aab = \ ~ 1}- For later calculations it turns out to be 
useful to know the number of elements in these sets. 

i^Aaa = 2a, h-1 e Aaa, 

i^Aab = 2a, h-1 4 ^ab, a <b, 

H^Aas = a, h-1 4 Aas, (19) 

#Ass = i^As's'=n/2, h-1 e Ass,As's', 

i^Ass' = {n-2)/2, h-1 i A,,,. 

We assumed a^h ^{s, s'}. For odd n these results have to be slightly modified. 

The fusion angles can be divided into three classes. The first one involves the particles s 
and s' . 
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= ul, = h-2a, ul, = ul,, = h/2 + a. (20) 

The second and third class do not incorporate s and s' and are distinguished by a relation 
of the three particles in the the fusion process. We have 

u'^fj = h — c, = h — b, ul^ = h — a, if a + b + c = h, 

(21) 

ulfj = h — c, ul^ = h — a, Uac = if a — b + c = 0. 

Up to now we had to review a lot of known things and to introduce a bulk of notation. 
We are now in a position to come to the main part of this paper. 

There might be several ways to write down the part of the full form factor which contains 
poles and zeros in the physical region. We prefer a slightly modified version of what has 
been presented in |]12[ . This version is based on a factorization of the pole part of K while 
the part containing zeroes is the quantity to be determined by the form factor equations. 
We take Xj = e^' and define the follwing vectors 

x^'^) = (xf \ , ke{l,2,...,n-2,s, s'}. (22) 

The Nk denote the number of particles of type k in the form factor (cf. (^)). We can then 
write down the following parametrization: 

(23) 

The factors l/{x[''^ + x^-^'') parametrize the kinematical singularities. The product in 
front of them is already chosen for D!^^ with n even, i.e. all particles are self conjugate. 
Obviously we can write down K without any difficulties for odd n as well. In what follows 
we give the explicit formulas for the cases with even n only. 

The objects Wki{xf'\xj'^) do parametrize the fusion poles. In [T^ it was shown that a 
particular ansatz for them does lead from (|^) and (0) to a polynomial equation for the 
Q's at least for A^^^-ATFT. We take this ansatz in a slightly modified form. 

Wkiix\''\xf) = n (4'^ - f^'^+'xf ) (xf ) - n-'-'xf), n = e*-/^ (24) 
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In the ansatz ( p3D above Q[Ni,N2,...,n^,] is assumed to be a polynomial which is due to the 
structure of ATFT symmetric at least in each of the vector components of any x^''^. Of 
course this assumption has to be justified case by case especially in the presence of higher 
order poles in ATFT. Since our form factors are required to be Lorentz invariant it is 
straightforward to calculate the degree of the polynomials Q. A Lorentz transformation 
consists simply of a linear shift of the rapidities. The denominator of (p3D is not Lorentz 
invariant. However its degree determines the total degree of Q to be: 



deg QiN„...,N,] = YTuZl (iVfc(iVfe - mk - I) + 2Nk{N, + Ns^)k) + 4E^=? Er=t+i kNM 
+ ^Y.k=s,s'Nk{Nk-l) + NsNs\n-2). 

(25) 

The partial degree, i.e. the maximal degree of a variable x\ in Q can be calculated only 
after having established the recursion equations for the polynomials below. In general 
this leads to a quite complicated formula. We treat the problem of calculating the partial 
degree for the special case of D^l"* in section 5. 

It is now clear that we have chosen an ansatz for the form factor in such a way that we 
have to solve only for the polynomials Q. However, as can be read off from (^) these 
objects do have a quite high degree already in cases when only a few particles are present. 
This makes the general solution technically quite difficult. 

We can now plug the above ansatz into the residue equations to obtain recursion relations 
for the polynomials Q. In doing this we first make the following observation 



X Qa, [^„...,^,](X(«),XW,X«, . . .X(^')) g^(xW, . . . ,x(^')). 

L 1 ! ■ ■ ■ ) \ 

(26) 

This equation indicates how the two particles a and h undergoing either a kinematical or 
a bound state process are separated from all other particles in the form factor. The last 
factor in the product does not contain any dependence on the particles a and 6, i.e. the 
coordinates x^"''^ and x*^^-* do not appear in that expression. A similar expression can of 
course be written down if we want to seperate only one particle. 
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By using the structure of the index sets Aab of the exact S-matrices for D^^-'-ATFT and 
introducing the following piece of notation to make expressions a bit more transparent 

[r]f ) = X - fi'- xf \ A; G {1, 2, . . . , n - 2, s, s'}, ^ = 1, 2, . . . iV^, (27) 

we arrive from the kinematical residue equation (P) at the following recursion equation 
which links an + 2-particle polynomial to an A^-particle polynomial. 



gaa[^....^.,](-a:('^\a;W,x«,x(2), . . . ,x(^')) = -,(-1)M.(xM)4«-i 1Ua^^^(^.. ) 

X (nti Ul\ n.a.. [-r-l + [r + it UreA^, [1 - r - [B - r - 1^ 

- nti n£=i UreA^, [-h + r + Ip l-r - l]f ^ UreA^, [r - 1 + B^ [r + 1 - B^) 
X Q[;v....;v.,](x«,x('\...,x(^')). 

(28) 

Obtaining a closed recursion relations for the polynomials Q coming from the first order 
bound state equation ([T^) is not so straightforward. It was mentioned above that the 
fusings in D^^^-ATFT can be devided into three classes (^O]), (^if) . In the second class of 
(pl|) the three particles entering appear to be unsymmetric in the sense a — b + c = 0. 
In order to calculate the desired recursion relation these three classes have to be treated 
seperately and especially the third class carefully. We also have to take into account the 
number of elements in the sets Aab (|19D. The first order equation for the process a + b ^ c 
is then 



X nti Ul\{Ur.A^. [1 + r + ulf\-l + r + ulf\-ul + r-B + lf\~ul + r + B- l]f 
X Ur.A,, [-< - r + lf\-ul - r - lf\-r + ul + B- lf\-r + ut,-B + ifP) 

he 

X (5c[7Vi...iv,,](a;,x(^\...,x("')). 

(29) 
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The last two lines should be compared with the explicit expression for A given in ([T3|). 
u'^fj etc. are the fusion angles for the possible processes which involve the particles a, b, 
and c. We used the abbreviation 9^^ = + mJJ^, — h. For self-conjugate theories it holds 
of course that G^f, = u^^^. 

Even though we do have negative powers in the above expression the equation is a poly- 
nomial recursion equation. It can be checked that these negative powers are canceled 
for any ATFT. We have, however, chosen to write them in order to make the remaining 
products a bit more handy. 

We note already at this stage that the recursion coefficients in ( P^D and (^) can be 
expressed as a sum of elementary symmetric polynomials. We will treat this point more 
explicitly in section 5 when we calculate some solutions to these equations for the d'^^ 
case. 

In order to be able to obtain solutions to these equations without going into tedious case 
by case calculations it would be desirable to have a Lie algebraic interpretation of the 
two recursion equations for instance in the sense of [0. We were informed that an 



equation similar to ( pUj ) can also be obtained for the Q-polynomials arising in A^^-ioim 
factors. 

It is important to note that (^) is written if the particles are in order. This means 
that we impose an ordering in the set of particles {1, 2, . . . n — 2, s, s'} by 1 < 2, ... 1 < 



s', . . . , s < s'. Then in (|29| ) we have a < b. 

If we change this order of the particles in the process, i.e. we consider b + a —>■ c rather 
than a + b ^ c the coefficient in (p9D will change. This change is simply obtained by 
exchanging the first two arguments of Q and replacing Q in the coefficient. Since 

both processes are physically equivalent this artifact leads to an identity for one particular 
polynomial Q at specific values of the arguments. This will become more transparent in 
section 5. 

Let us now come to a discussion of the form factor equations in the presence of higher 
order poles of the S-matrix. The content of equation (^9]) can be stated graphically by 
the following diagram which should not strictly be intrepreted in the perturbative sense. 



b 

FIG, 1: First order pole diagram in the form factor 



The bubble in this figure depicts the local operator in the conventions of 
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Let us consider the third order foward channel poles of the S-matrix first. These poles 
correspond to fusings which are already covered by the fusion angle analysis at the be- 
ginning of this section. For our purposes this means the following. The perturbative 
analysis shows that in general several Feynman diagrams are needed in order to explain 
the third order poles of the S-matrix. But the S-matrix in our case is known to satisfy 
the bootstrap equations. Hence, for the forward channel fusings at the position of a third 
order pole of the S-matrix we can in analogy to fig.l draw the following effective picture. 



D- 



FIG, 2: Effective picture for the third order pole 

The circle at the position of the fusing indicates that we do not need to know about the 
perturbative processes which lead to the production of the physical particle running into 
the local operator. 

This means that we should be able to derive equations for the Q-polynomials correspond- 
ing to fusings of the type in fig. 2 in analogy to what has been done for the first order pole 
fusings. 

A detailed analysis shows that this is can consistenly be achieved. The result for a fusing 
a -|- 6 ^ c at a third order pole of the S-matrix is exactly equation ( pQ]) where we only 
have to replace the residue appearing in the denominator of the right hand side of (|29|) 
by the coefficient of the leading order singularity of Wab{x' , x)~^ . 

In other words, due to the bootstrap properties of the S-matrix the fusings corresponding 
to third order poles of the S-matrix can be treated on exactly the same footing like 
the first order ones and the corresponding recursion equations for the Q-polynomials are 
structurally identical! 

We conjecture that the same phenomenon is true for the higher order forward channel 
fusings in the ATFTs coming from the E-series as well. 

The problem of deriving equations for the form factor in presence of higher odd order 
poles of the S-matrix has already been addressed in the context of Ising models in |TT|. 



The question arises what happens in the case of the second and fourth order poles of 
the S-matrix. For a discussion of this question see also |]TT|. These poles are known 
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T5[] not to participate in the S-matrix bootstrap. It is therefore not possible to draw 
effective diagrams hke the one in fig. 2 for these cases. One might therefore be tempted 
to think that these poles do not lead to any significant contribution on the form factor 
level. However, this is not true. We will see in the next section that it is possible to 
derive a consistent equation for the Q-polynomials corresponding to the second order 
pole processes of the S-matrix in D4^''-ATFT. Moreover, it will turn out in section 5 that 
these equations are needed in order to constrain the solution spaces of the Q-polynomials. 
Hence it seems that for the even order poles the processes which are needed in order 
to explain the corresponding S-matrix element have an influence on the solution spaces 
of the form factors. It is, however, not clear how to derive in general equations for the 
Q-polynomials in these cases even though it might in principle be clear how the equations 
for the full form factors can be formulated (see e.g. |]1T], 0, ^ and (pQ])). The main 
reason for this mainly the fact that in D^^^-ATFT (in contrast to the A^^ cases) we do 
not know a priori how many Feynman diagrams do contribute to an even order pole of 
the S-matrix and what kinds of particles are involved in a particular diagram. Therefore, 
up to our present knowledge the contributions of these diagrams have to be studied case 
by case. As mentioned above we are going to treat the case of the second order pole in 
D^^ in detail in the following sections. 



4 The case 



We are now going to discuss the form factors in the simplest D^^-ATFT. As mentioned 
above the /^^^''-ATFT describes four particles which we label 1, 2, s, and s'. The mass 
of particle 2 is rrih = ^/^m while the masses of the other particles are the same with 
rrii = \p2m. This means that we can devide the particle content into light particles 
/ G s'} and one heavy particle h = 2. This kind of symmetry follows from the 

sjTiimetry of the Dynkin diagram of D4 (see e.g. fl^). In this section we are going to 
present the form factor equations of this particular model. We will see that the symmetries 
of the Dynkin diagram are maintained at the level of the form factors. 
Even though this model is quite simple compared to the general D^^^-models it already 
shows the generic features which are present in D^^^-theories. Referring to the S-matrix 
we do have a third order fusion pole for the process h + h ^ h and a second order pole 
for the process I + h —* I' + I" . However, the S-matrix does not have fourth order poles 
which makes life a bit easier. 
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We are going to apply the machinery presented in the previous section. First we use ( ]25|) 
and obtain for the total degree of our Q-polynomials: 



degQ[N,M,N,] = |{N^{N^-l) + Ns{N,-l) + NANs'-l)) + lN2{N,-l) 
+ 2 {N,Ns + N,Ns' + NsN,,) + AN2 (iVi + iV, + iV,,) • 

(30) 

Using the structure of the S-matrices we obtain the following equations for the kinematical 
poles. 

n , J -y T y(1) - „3 -3i (~1)^' 

in^M?[^f[-B]?[B - 2f[-B - Af[B - 6f 

nfii([3]f )[5]f [-5 - Mf\B - 3]f [-5 - 3]f\B - 5]f ) (31) 

n.=.,z" nltMXn-B - 2mB - 4]L') 

M-[-r])g[^,,,v„7v.,iv,,](x«,...,x(^')). 

The notation in the last line indicates that on the right side of the minus sign appears a 
polynomial structurally identical the the one on the left side with all the entries in the 
brackets replaced by their negative values. 

In this expression / can be either 1, s, or s' and therefore shows the symmetry of the 
theory on the form factor equation level. 
For the kinematical pole of particle h we find 

^hh[Ni,N2,Ns,N^,][—X,X,X.^',...,X.'^ ') — X — pSnTT-Tj— 

hh ^ ^ 

(n.=M',P n£i[l]f ^([3]f ))^[5]f [-i? - lf\B - 3f\-B - 3f\B - 5]f ^ 
nfii([2]f V([4]f Vh^lf ([-5 - 2]f^f[-B - A]f\B - 2]f\[B - 4]f - 6]f 



-r]) g[iVi,7V2,7V.,iV,,](xW,.--,x(^')). 



(1) 

(32) 
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To get the recursion equations for the first order bound-state residue equations we intro- 
duce the following notation to make the formulae as transparent as possible. 

{r}f ) = [r + lf\r - lf\-h + r + 1 - Bf\-h + r - 1 + Bf\ (33) 

Here h stands for the Coxeter number which is equal to six in the present case. 

We can derive the following recursion relations in D^^^ where we put some constant factors 
appearing in (^) into a constant H for the particular process. 
First we get for I + I ^ h: 



QUIN,...N,] {XQ, x(l), X(2), x(^), x(^')) = 

Ha U^,^^Ul\[Qf^ g,[^,...^,](x,x«,x(2),x(^),x(«')) 



(34) 



The next processes to be considered are I + h ^ I and I + I' ^ I" . Using (|29|) we arrive 
at the following recursion relations 

X llj=l J^jll) lli=l ^p) lli=llC/i lli=llC>/i ^l[Ni...N^,][X,X' ' , ■ ■ ■ '), 

(35) 

and 

Qu'[N,...NA^n',xn-^^^'\ . . . ,x(«')) = Hw x' 

(36) 

X nfj, ^ n£ u^mf" nS mPQnN....N^, . . . , x(^')), 

respectively. The next process which occurs is h + h —>■ h and is of third order with respect 
to the S-matrix which means that in the ansatz for K (ESf) a pole of second order appears. 



We evaluate according to the rules for third order forward channel fusion poles and arrive 
at the following recursion relation 

QM[7Vl...iv,,](a;^^^a;^^~^x(l),...,x("')) = Hhhx'^ 

X lli=i j4j(h)jg](/i) PJi [lk=l,l',l"i'^Ii V/i[Afi...Af,,]l2;,X^ ^ . . . ,X^ 
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One recognizes that structurally there is no formal difference between the recursion re- 
lations which arise from a first order pole and the one coming from a third order pole 
respectively. 

Let us now address the problem of the second order pole which occurs aX 9 = m/2. 
The diagram corresponding to this pole and the corresponding cut-diagram which leads 
to a form factor equation is as follows (see also |]lT], Q). 




FIG, 3: Diagram for tlie second order pole 



Our goal is to derive a relation for the Q-polynomials arising from the second order 
process. We take the full form factor (^) together with our parametrization of the pole 
part (0). Of course this ansatz can be taken for granted only after having established its 
consistency with the second order equation. This means that this ansatz is sufficient in 
the presence of a second order pole if and only if we derive a polynomial recursion relation 
for the Q's. This means we do have to verify that the Q's are actually polynomials. 
First we have to calculate the momenta of the particles in fig.3. These are at the value of 
the pole 9i = d + m/2, Oh = 6, and 6i, = Ov> = 9 + m/Q. 

Using the explicit form of the minimal form factors [O, 11^ we can deduce the following 



identitiy in D 



(1) 



^' pmin 

nn 



6f^+m/Q)F{^f{d-er+m/Q) 



(fc) 



k=l i=l 



Fin 



n n , .(h) ' 

i=i '=1 \^/+{e^) 



) + m/2) - 0. 

If we then write the equation for the full form factors in the following way 



(3^ 



-Z reS6)'=6l+i7r/2 ^lhdi...dn 



r, 6,01,..., On) = rf.i„rU.F,i"d^...dA0+^^/Q, e+nx/Q, Oi,..., ft„), 

(39) 
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which is in accordance with the two particle form factor equations in an Ising model [[TJ 
and with [21|, we arrive at another polynomial equation for the Q's 



X nfi, ^ Y^A ^ n£[6]f nS [6]?"^ Qv.,n.,...nA-^ x(^), . . . , x^'). 



(40) 



We can verify the following identity, which should be compared with an expression which 
occured in |jl2| in connection with the second order poles in A^^"* 



This identity is necessary in order to guarantee the consistency of the solutions of P0| ) 
with the other equations in D^l^ due to the presence of minimal form factors at particular 
values in the equations (^) and (|29|) . 

It is remarkable to note that even (^) can more or less be treated on the same footing 
like the first order equations for the form factor by realizing that if the fictitious vertex 
(cf. d^)) Vih = —i''Tes''e=in/2Sih{0) is equal to the expression Ti„i„T\,i„ which occurs in 
(^). The quotation marks at the residue mean that we have to evaluate at the leading 
order of the singularity. 

Moreover (1^) looks very similar to the ones coming from first order processes. In par- 
ticular one might have noticed above that the particles running into the local operator 
acquire a factor [6] = [h] in the recursion coefficient. This also happens here where the 
two light particles running into the local opeartor just have this factor in the recursion 
coefficient. 

To all the bound state equations the remark on the order of the particles which has been 
made after stating the general equation (pOf ) applies. This means that we get additional 
relations for one and the same Q-polynomial at specific values of the arguments ! 

In the next section we will give some details on the calculation on solutions to the recursion 
relations derived above. It turns out that we do neccessarily need equation (|40|) which 
originates from the second order pole of the S-matrix. It is needed to reduce the degrees 
of freedom of the solutions. 



This is an important fact because it is known and it was mentioned above that the 
second order diagrams do not really contribute to the S-matrix bootstrap but they are 
needed for the corresponding form factor bootstrap ! 
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Let us briefly comment on second order poles in general. In principle we can expect an 
equation for the entire form factor of the kind (^). Even though the position of the second 
order poles of the S-matrices in the D-series are known there is no closed expression in 
the literature for the particles actually participating in the second order process. We have 
checked for the Dq case that an equation of the type ( pUj ) gives a result comparable to 
(0) for second order processes. However, since ( ^OD has a quite nice structure one might 
expect to be able to write down a closed recursion formula for all admissible second order 
processes in the D-series. 



5 Some notes on solutions in the case 



Now we are coming to the construction of solutions to the recursion relations derived in 
the previous section. 

It was already pointed out that due to the ansatz for the singular part K of the form factor 
(^) the polynomials Q[Ni,N2,Ns,n^,]{'x^^\^^'^\^^'^\^^^ ^) are symmetric at least in each of 
the components of one particular coordinate vector with k = 1,2, s,s'. It might 
therefore be useful to expand Q in symmetric polynomials. A basis in the space of these 
polynomials which is commonly used in form factor calculations (see e.g. §|, 01) are the 
so called elementary symmetric polynomials e^"^ However, even though we are going to 
use them in the calculations in this sections we have the impression that they might not 
be the appropriate basis for the problem because the actual calculations become quite 
tedious when using this basis. We will comment on that problem in the final section. 



Let us introduce some simple facts about elementary symmetric polynomials [Q. As 
functions of n variables xi,X2, ■ ■ ■ ,Xn they are defined by the following expression 

n n 

n(l + tx,) = 5:eWr. (42) 

i=l r=0 

Here the superscript on the e's refers to the number of arguments. Explicitly the e's are 
given by Cq = 1, e^"^ = xi + X2 + ■ ■ ■ Xn, . . ., e^^^ = xiX2 ■ ■ ■ x„. 

Comparing ( pTj) with (^2]) one recognizes at once that all the recursion coefficients in the 
equations for the Q-polynomials can in principle be expressed in terms of the elementary 
symmetric polynomials. 
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The next entity we need is the concept of a partition A of a positive integer m. For our 
purposes this is a finite sequence of non-negative integers arranged in decreasing order 
A = (Ai, A2, . . . , K) with Ai > A2 > . . . > Ar- Characteristic quantities for A are its length 
( /(A) = r in our example ) and its weight |A| = Aj (which is equal to m in our 
example) . 

Having these concepts we define an elemtary symmetric function associated to a partition 
A by 

= 4^--ei^ (43) 

It is then clear that the weight of the partition gives the total degree of E^^^ while its 
length gives the maximal power of a particular variable Xi in The latter is called 

partial degree. 

Since we have recognized that the coefficients in the recursion equations are in fact ele- 
mentary symmetric polynomials we can read off these equations in how much the degree 
of a variable x-^^ is changed when performing one step in the iteration of the polynomials 
Q. Since the recursion equations are linked it is not quite a straightforward task to cal- 
culate the partial degree of the coordinates of {x^^\x^'^\x.^'^\x^^'^). Given the polynomial 
Q[Ar(,Ar,,,Ar;„,Ar,j(x'^'\ x'^''\ x'^'"\ x*^^^) where Z, /" G {1, s, s'} we find for the maximal length 
of the partitions A^'^ corresponding to particle I 

/(A(')) < 3{Ni - 1) + 2{Ni, + Ni.) + 4Nh. (44) 

The partial degree of I' and I" resp. can be obtained by just replacing the letters in (Q). 
For the heavy particle we get 

/(AW) < 7{Nf, - 1) + 4{Ni + N, + N,,). (45) 

This result refiects again the symmetry of the Dynkin diagram at the level of form factors. 

Due to the presence of four different particles in the theory we have to associate an ele- 
mentary symmetric polynomial to each of the particles seperately. We therefore introduce 
the following notation 

A = (A«|A(2)|A(^)|A(^')). (46) 

There are many different ways to choose initial conditions for the symmetric polynomials 
Q, see e.g. 0. From the degree-formulas (^), ( |30D one can see that the Q-polynomials 
with only one argument are of zero degree which means that they have to be constant. 
We therefore choose the following initial values for the Q's 
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Q[1,0,0,0] = ^(1), Q[0,1,0,0] = ^(2), Q[0,0,1,0] = ^(s), Q[0,0,0,1] = ^(s')- (^7) 

With this data we can compute the two-particle polynomials. With d'^''^''^'^ we denote 
constants which can not be determined by the recursion equations. 
If two light particles are present we find 

Q[2,o,o,o](^i'\ 4'^) = i^^ii^m - 3cf''°'°'°l)£;(2i|o|0|o) + cP'°'°'°l£;(i3|o|o|o). (48) 

We have written down the polynomial for two particles of type 1. The polynomials for 
particles s and s' are structurally identical. One has only to change the numbering at Q 
and the £"s etc. to the appropriate place. 

The next polynomial corresponds to two different light particles. We again write only one 
member of this family, the other ones can be obtained by symmetry. 



<5[i,o,i,o](4'\ 4'^) = c^'''''''\E^i^mo) + ^(o|o|i^|o)) + {Hu>ns> + c^'''''''^)E^^\om- (49) 

The polynomial for two heavy particles is already of degree 7 and it is on the two particle 
level not possible to eliminate most of the constants. 

I (2) (2)\ [2,0,0,0] rn , [2,0,0,0] j-, 

V [0,2,0,0] l^^i ) — Ci £/(o|i7|o|o) + C2 -c^(0|2l5|0|0) 

, [2,0,0,0] IT, , [2,0,0,0] IT, 

+C3 -C/(0|22l3|0|0) + C4 -C/(0|23l|0|0), ^5Q^ 

[2,0,0,0] [2,0,0,0] [2,0,0,0] [2,0,0,0] _ rr 

The case of one light and one heavy particle is interesting because we can see there that 
the equation arising from the second order pole of the S- matrix is needed in order to reduce 
the number of degrees of freedom of the solutions. The basis for this particular case which 
is a polynomial of degree 4 is spanned by E(i4|o|o|o), £'(i3|i|o|o), £^(i2|i2|o|o), £^(i|i3|o|o), and 
-S(o|i4|o|o)- If we label the corresponding constants consecutively for these polynomials 
we see that the first order equations reduce to three degrees of freedom (one of them is 
actually reduced by the symmetry properties of the Q-polynomials at particular arguments 
mentioned above), while the "second order equations" remove another two degrees of 
freedom. We then find 
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^[1,1,0,0] ^ l(-4MA0]^p(^^2^]-3)), 4^'^'°'°] = (-i-^)c?'^'°'°l-^(i/-P), 

[1,1,0,0] _ / 1 _ n-:^Jl, 1,0,0] _ Q-3 1^3 p [1,1,0,0] _ _i [1,1,0,0] 

(51) 

where we have set H = HihTii^i) and P = sS'^'^'^'^^ + Huirii^s')- 

It is clear already at this point that the structure of the solutions becomes quite compli- 
cated already at the two particle level due to the high polynomial degrees. 
Of course it is straightforward, however tedious, to compute the polynomials for higher 
particle content as well. For example the polynomial (5[2,i.o,o] contains almost thirty 
elementary symmetric polynomials. We do prefer not to list the higher cases in this 
paper. 

We would like to add the following remark. In (^4D and (|^) we already obtained a 
result on the partial degree of the Q-polynomials. It can be taken for granted that all the 
partitions for a specific Q which are constrained by the partial degree do contribute to the 
Q-polynomial. This means that we will not find a case where the coefficients c[.'^'^''=''^l for 
the corresponding basis function do vanish. In other words all the elementary symmetric 
basis functions obtained by this general rule do and must nonvanishingly contribute to a 
consistent solution Q[a,fe,c,d]- 

We do believe that a basis using elementary symmetric polynomials is not the appropriate 
one for explicitly calculating the Q-part of the form factors and therefore would like to 
postpone a detailed discussion of the higher order polynomials. We are going to comment 
a bit more on that point in the next section. 



6 Discussion 



By using the form factor bootstrap and a specific (but generic) ansatz (^) for the singular 
part of the form factors we have been able to derive closed recursion relations for Q- 
polynomials within this ansatz. One of these equations arises from the kinematical residue 
equation of the form factor bootstrap. The other one which is at first a consequence of the 
first order bound state residue equation of the form factor bootstrap turned out with some 
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slight modifications to be applicable for the fusings in D^^-'-ATFT which come from third 
order poles of the S-matrices as well which is a consequence of the S-matrix bootstrap. 
We did comment on the consequences of the second and fourth order poles of the S- 
matrix on the Q-polynomials. Since at present there is, in contrast to the A^^ cases, 
no rule saying which particles do participate in the Feynman diagrams contributing to 
the even order poles of the S-matrix we cannot give a general closed equation for the 
Q-polynomials in these cases. However, it is clear that at least some of them are needed 
in a particular theory in order to reduce the number of degrees of freedom of the Q's. 
We treated D^^-ATFT special case of our general analysis. We have shown how 
the symmetries of the Dynkin diagram are in a sense maintained in the equations for 
the Q-polynomials in this case. Moreover we treated the influence of the second order 
poles of the S-matrix on the form factors in detail and derived another equation for the 
Q-polynomials. 

In attempting to calculate solutions in the D^^ case we realized that due to the high 
polynomial degree of Q the solution spaces are very complicated when a basis for Q 
consisting of elementary symmetric polynomials is used. 

In order to find solutions without going through tedious computations it would first be 
interesting to find a Lie algebraic interpretation of the recursion relations. Second one 
has to think of another polynomial basis in the solution spaces. A direct guess would 
be to use Macdonald symmetric polynomials [^. This is because these polynomials are 
symmetric polynomials which depend on two parameters. This fits well into the affine 
Toda case since we do have exactly two parameters in the recursion equations. One of 
them is f2 = e*'^/'^ and for the other one we can take ^l^ with B being the effective 
coupling (^. Due to the weak-strong duality B ^ 2 — B |15] and the symmetry of the 



Q-polynomials which was pointed out in the discussion after equation we have some 
hint that the Macdonald polynomials might be a proper candidate (see ||2^). Also their 
connection to root systems of Lie algebra (see e.g. |2^) might be of use for the present 
problem. 

Let us make one last comment on the deformation aspect of the solution spaces. In 
form factors of the scaling Lee- Yang model were calculated. If one considers the closed 
solution of PI for the Q-polynomials one finds that a solution to the n particle case is 
exactly one skew Schur function 



Qn ~ sx/^{xi, . . .,Xn), A' = (2n - 2, 2n - 3, . . . , n - 2) , /i' = (2n - 4, 2n - 6, . . . , 0), 

(52) 

where A', /x' are the conjugate partitions. 
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The scaling Lee- Yang model can in some sense be considered as a special case of the 
sinh-Gordon model (more or less by setting B ^ 0). Or the latter one is a deformation of 
the first. In it was shown that for the sinh-Gordon form factor problem it is possible 
to find a closed solution for the Q-polynomials as well. The result is 

Qnik) = detM.jik), M,,{k) = [i-j + k]eflj, (53) 

where M is a (n — 1) x (n — 1) matrix and [n] = sin(ni?/2)/sin(i?/2). This solution is a 
then a certain deformation of the skew Schur function above. 

It would be interesting in the case of D^^^ theories if in certain limits models exist for 
which a closed solution for the Q's can be constructed easily which then would allow for 
hints for the generic cases. 
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